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Abstract 

We demonstrate that a new type of analysis in heavy-ion coUisions, based on an event- 
by-event analysis of the transverse momentum distribution, allows us to obtain information 
on secondary interactions and collective behaviour that is not available from the inclusive 
spectra. Using a random walk model as a simple phenomenological description of initial state 
scattering in collisions with heavy nuclei, we show that the event-by-event measurement 
allows a quantitative determination of this effect, well within the resolution achievable with 
the new generation of large acceptance hadron spectrometers. The preliminary data of the 
NA49 collaboration on transverse momentum fluctuations indicate qualitatively different 
behaviour than that obtained within the random walk model. The results are discussed in 
relation to the thermodynamic and hydrodynamic description of nuclear collisions. 

1 Introduction 

The ultimate goal of present heavy ion experiments at the CERN SPS and similar studies 
foreseen at future collider facilities is the production and characterisation of an extended 
volume of deconfined quarks and gluons, the quark gluon plasma (QGP). The QGP has 
been predicted as the equilibrium state of strongly interacting matter at sufhciently high 
temperature and density 0| . A variety of possible signatures for the transient existence of 
a deconfined phase of matter in the course of nucleus-nucleus (A+A) collisions has been 
proposed theoretically and studied experimentally (see for a recent review) . 
An important aspect of these studies is answering the question to which extent and at 
which stage of the interaction thermal and chemical equilibration is reached in nuclear 
collisions. Recent experimental results show a significant change of the shapes of particle 
spectra and relative particle yields when going from nucleon-nucleon (N+N) interactions to 
central collisions of heavy nuclei . These observations find a natural interpretation within 
statistical (equilibrium) models assuming thermal and chemical equilibration of produced 
matter (hadron gas 0, D or QGP Q). However there are numerous attempts to describe 
the data without invoking equilibration and assuming that the observed effects are only due 
to initial, non-equilibrium, scatterings. In the case of inclusive observables equilibrium and 



non-equilibrium approaches lead to a qualitatively similar behaviour and therefore they are 
difficult to distinguish, if the absence of the caluclable theory in the domain of soft QCD 
processes is taken into account. 

It was realized Q that a study of event-by-event fluctuations may serve as a tool to 
measure equilibartion in nuclear collisions and therefore allow to distinguish between equi- 
librium and non-equilibrium models. In this paper we follow this approach and consider 
transverse momentum fluctuations in non-equilibrium, initial state scattering models and 
in equilibrium, hydrodynamical models. 

It is experimentally established 1^ that the inverse slope of transverse momentum dis- 
tribution increases when going from p-fp interactions to central Pb+Pb collisions. The 
strongest effect is observed for heavy hadrons. The equilibrium models relate this obser- 
vation to the presence of transverse hydrodynamical expansion of the matter Q. The size 
of this effect increases with the size of the colliding nuclei and is most pronounced in the 
transverse momentum spectra of the heaviest particles. Such models, which apart from the 
choice of an appropriate flow profile, can be characterized by two free parameters, the tem- 
perature at freeze-out and a flow velocity parameter, have been shown to fit the measured 
inclusive spectra ^, ||] . 

However, it has been demonstrated that the same spectra can also be fitted by a different 
class of models, which do not require the existence of equilibration [|o[ |ll|. These models 
are triggered by the observation of a broadening of transverse momentum distributions al- 
ready in proton-nucleus (p-|-A) interactions ('Cronin-effect'), which is commonly attributed 
to multiple scattering of the incident nucleon and its remnants in the target nucleus Q . 
Several authors have parametrized the initial state scattering, which effectively translates 
part of the longitudinal momentum of the incident nucleon into transverse momentum, in 
terms of a 'random walk' in a transverse kinematical variable such as a transverse rapidity 
pO| or transverse velocity . The parameters of the random walk can be fixed by a fit 
to p-|-A data leading to a fit of A-l-A spectra that is essentially parameter free, if the 
inverse slopes from N-t-N interactions are used as starting values. Keeping this difference 
in mind, the initial state scattering or random walk models provide a description of the 
data that is of similar quality as that of the hydrodynamical models. A study of inclusive 
particle spectra alone is therefore clearly not enough to understand whether the produced 
particles originated from a source having a sufficient amount of secondary interactions to 
develop collective behaviour such as hydrodynamical transverse flow. These two approaches 
can however be distinguished by analysis of two pion correlations [p^ and event-by-event 
fluctuations as shown in this paper. 

The paper is organized as follows. In Section 2 we recall the method introduced in Ref. 
and used further in the paper to study event-by-event fluctuations. The fluctuations in 
an initial state scattering model are calculated in Section 3. The discussion of fluctuations 
in thermal models is given in Section 4. Summary and conclusions close the paper. 



2 The Measure of Event-by-Event Fluctuations 

Event-by-event fluctuations in nuclear collisions are usually dominated by the trivial vari- 
ation in impact parameter from event to event and the purely statistical variation of the 
measured quantities. A statistical method that allows us to remove these trivial contribu- 
tions and to determine the dynamical event-by-event fluctuations has been proposed in m . 
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Following this reference we define for every particle i: 

Zi = Pt, - Pt, 

where is the mean transverse momentum of accepted particles averaged over all events 
(the inclusive mean). Using Zi we calculate for every event 

N 

1=1 

where N is the number of accepted particles in the event. With this definition we obtain 
the following measure to characterise the degree of fluctuation in the pr distribution from 
event to event: 

where (N) and {Z'^) are averages over all events and the second term in the r.h.s. is the 
square root of the second moment of the inclusive transverse momentum distribution. 
The physical motivation for studying was given in experimental data on N+N in- 
teractions show that particles in these collisions are not produced independently One 
observes large scale correlations that lead to, e.g., a correlation between the event multiplic- 
ity and the average pr of the particles. This property can be used to probe the dynamics of 
nuclear collisions by measuring to which degree the correlations present in N+N interactions 
are changed when going to p+A and A+A collisions. 

For this purpose, as a measure of fluctuations has two important properties. For a 
large system (i.e. a A+A collision) that is a superposition of many independent elementary 
systems (i.e. N+N interactions), has a constant value that is identical to that of the 
elementary system. In other words if the central Pb+Pb collisions were a incoherent super- 
position of N+N interactions, the value of would remain constant, independent of the 
number of superimposed elementary interactions in a single event and its distribution in the 
studied sample of the events. If on the other hand the large system consists of particles that 
have been emitted independently, assumes a value of zero. Thus ^pj. provides us with a 
scale characterising the fluctuations in nuclear collisions relative to elementary interactions 
at the same energy. 

One should expect that $py is sensitive to the initial state scattering effects considered in 
||lO|| and through additional fluctuations in the transverse momenta of effective N+N 
sources superimposed in these models. As the amount of rescattering is directly related to 
the size of the colliding nuclei, here one expects a smooth dependence of on the nuclear 
atomic number. Finally, if the underlying physical mechanism for A+A collisions corre- 
sponds to some different type of correlations between the (average) transverse momentum 
and multiplicity such as discussed, e.g., in the context of hydrodynamical models (l4[ |l6) 
and further in , then one could expect an abrupt change in $py , when crossing a bound- 
ary in the size of colliding nuclei and/or collision energy, where the new physical mechanisms 
beyond those in elementary interactions are already at work. 
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3 Event— by— Event Fluctuations 

in Initial State Scattering Models 

The motivation for applying this method to test the initial state scattering models is that par- 
ticles in these models are not produced independently. They are emitted from sources that, 
after an initial scattering, are moving in transverse direction. This transverse source move- 
ment is superimposed on the particle spectrum emitted by the individual particle source, 
where apart from the overall boost in transverse rapidity each N-|-N-like source is character- 
ized by the same correlation of transverse momentum vs multiplicity as for N-l-N interactions. 
Thus the pattern of particle emission changes. A new source of fluctuations is added (ini- 
tial state scattering) and we expect that the fluctuation measure from Eq. (§) will be 
sensitive to its presence. 

To get an estimate of the size of this effect we use as reference a simple parametrization 
of p-l-p data at ^/s « 20 GeV/c as proposed in Ref. 0. We restrict ourselves to a study 
of negativly charged hadrons, which for this energy are mostly (> 90%) pions. The trans- 
verse momentum, px, probability distribution for events with negatively charged hadron 
multiplicity n is parametrized as: 

/inbr) = C„pTe-™-/^("\ (2) 

where C„ is a normalization factor, mx is a transverse mass (m^ = \/ {p^ + rn^), where 
m is the pion mass) and T{n) is the inverse slope parameter. The T{n) dependence on 
multiplicity is taken to be: 

AT 

r(n) = To + (n-l) — , (3) 

where parameters Tq and AT/ An are taken to be 173 MeV and -6 MeV in accordance to 
p-|-p data at ^/s w 20 GeV/c jEsj. The multiplicity distribution is calculated using the 
parametrization given in Ref. |19(| . 

We use a Gaussian rapidity distribution to describe the longitudinal momentum distri- 
bution of the produced hadrons. To check the sensitivity of our results, we use two extreme 
values for the width of the rapidity distribution, ay — 1.6 which corresponds to the value 
observed in p+p, and ay — 0.8, which is closer to isotropic particle emission in the source 
rest frame. 

Starting from this type of elementary particle sources, we then describe a collision of 
heavier systems following the model presented in [|lOj . In this initial state scattering model 
one takes into account the rotation of the collisions axis of projectile nucleons that have 
already undergone a collision with a target nucleon. The effect of multiple collisions is 
parametrized by giving the sources of produced particles ('clusters' or 'fireballs') a Gaussian 
distribution in transverse rapidity p |lO| : 

fABip) = [^]i/2e(-pV5iB), (4) 

using Sab = {Na + Nb — 2)(5^. Here Na and Nb are the numbers of interacting nucleons 
from target and projectile nuclei and 6 is the average shift in transverse rapidity per collision. 

As has been shown in , this model is able to quantatively describe the pT^-broadening 
observed in collisions of heavy nuclei, using a value for i5 of 0.22 that was fixed by fitting 
spectra from p+W collisions. 
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The result of our simulation based on the model described above are summarized in Fig. 1. 
It shows the evolution of with the width of the source transverse rapidity distribution 
(i.e. increasing system size). Only particles with | ycMS |< 3 and px < S GeV/c were 
used for the calculation of <f>py. The solid line in Fig. 1 indicates the reference value of 
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Figure 1: Dynamical px fluctuation measure for the random walk model with ay ~ 1-6 
and (Jj, = 0.8 as a function of 6ab- The soild line indicates the results without initial state 
scattering (a linear superposition of p+p interactions). 

^pj, « 4.5 MeV/c, which we obtain using our parametrization of p+p data. This value is 
the result of the correlation between average transverse momentum and mean multiplicity 
in the parametrization. 

For both values of ay we observe a strong enhancement of $pj, above this reference value 
for the random walk model, reaching 10 MeV/c for ay = 0.8 and 50 MeV/c for ay = 1.6 at 
a value of Sab = 0.44, which corresponds to central Pb+Pb collisions. 

The conclusion we draw is that the introduction of the random walk source movement does 
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not tend to wash out the correlations present in p+p interactions. On the contrary the 
correlation introduced by the transverse source velocity leads to an amplification of large 
scale fluctuations in transverse momentum even by a factor of up to 10. 

The preliminary results of the NA49 Collaboration show a qualitatively different 
behaviour. The value obtained for central Pb+Pb collisions at 158 A-GeV is significantly 
smaller than the corresponding value estimated for p+p interactions. This interesting result 
indicates that the non-equilibrium initial state scattering models are unable to reproduce 
the measured magnitude of the event-by-event fluctuations in central collisions of heavy 
nuclei. Note that similar conclusion was reached very recently by the analysis of two pion 
correlations within the same class of non-equilibrium models fl^ . 



4 Event— by— Event Fluctuations 
in Thermodynamic Models 

At the opposite extreme of the initial state scattering models without equilibration of the 
produced particles are thermodynamic approaches assuming global equilibration, with the 
system possibly undergoing a subsequent collective expansion until freeze-out. 
Initial global equilibrium is one of the basic assumptions of the Landau picture of nuclear 



collisions |21 , where the matter is fully stopped and globally equilibrated before the onset 
of expansion. 

In the Bjorken picture p2| initial termalization occurs at some proper time tq, setting the 
stage for the subsequent expansion. The basic difference between the two pictures is that 
in the Bjorken one the initial conditions for hydrodynamical expansion are assumed to be 
invariant under Lorentz boost in longitudinal direction. 



Leaving a detailed analysis of the event-by-event fluctuation pattern in equilibrium 
models to future publications, we comment on some points related to the above discussion. 
We have seen that the variable was designed in such a way that it shows the presence 
of correlations between the transverse momentum and multiplicity at the event-by-event 
basis. The starting point was the presence of such correlations in N+N interactions. At 
low collision energies the observed negative correlation can be interpreted as being due to 
energy and momentum conservation. The crucial question is then whether thermalization 
and hydrodynamical expansion can wash out such correlations. 

Let us start with the simplest example of an ultrarelativistic gas in thermal equilibrium. 
Then the number of particles is determined by the condition of thermal equilibration itself 
and the externally defined volume and total energy of the gas. This leads to a correlation 
between entropy (multiplicity) and temperature (transverse momentum) at the level of av- 
erage values {S ~ T^), when the external conditions vary. It is important to note, that the 
averaging here can be done with respect to a single system by measuring the quantities in 
question for a sufficiently long time. Moreover, the theory of thermodynamical fluctuations 
]23t shows, that the fluctuations of entropy and temperature in a given subvolume of a gas 
are correlated: 

{AS AT) = T, 

which for massless particles means a positive correlation between the multiplicity and the 
transverse momentum. Therefore, in a thermalized ultrarelativisitc system the particles 
are in fact not emitted independently and some definite correlation between the transverse 
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momentum and multiplicity takes place. Thus even in the case of full thermalization fluc- 
tuations due to the non-independent particle emission in heavy ion collisions should be 
expected. 

The physical origin of the fluctuations in the framework of hydrodynamical approach is 
the same. Following the above discussion one may distinguish between two different sources 
of fluctuations in hydrodynamical models. First of all due to fluctuations in the thermalized 
energy and volume occupied by the matter at the early stage initial conditions vary from 
event to event. In Ref. it was proposed to describe these early stage global fluctuations 
also in terms of thermodynamical fluctuations, where the role of the heat bath is played 
by the initial energy of the colliding nuclei. These fluctuations are sensitive to the number 
of effective degrees of freedom at the early stage and therefore, if measurable, can serve 
as a probe of the properties of the early stage matter. Finally thermal fluctuations in the 
elements of freezing-out matter take place . Here effective degrees of freedom are hadrons 
and hadronic resonances. As suggested in Ref. the magnitude of these fluctuations can 
be used to establish freeze-out conditions of hadronic gas. 

The exact form of the finally observed fluctuations depends on the type of expansion. 
Despite the considerable effort in investigating this problem [^4[ |l5[ |l^ further studies 
are needed to compare the predictions with the experimental data. 

We conclude, that in passing from the regime of incoherent N+N interactions to that of 
formation of thermalized matter and its subsequent hydrodynamical expansion, the corre- 
lation between the multiplicity and transverse momentum changes its origin. 
Thus we can expect changes in the value of the variable ^t- In the thermalized events 
showing a collective behaviour one has, roughly speaking, only one 'fluctuating unit', the 
whole system, whereas in the Cronin rescattering case the event is superimposed from a 
number of sources with fluctuating positions in the phase space. A reasonable guess is then 
that the relative "collective" fluctuations measured by $t should be smaller than in the 
Cronin rescattering case and the value of should thus drop, but further quantitative 
studies are clearly needed to settle this question. Let us also note, that in order to achieve 
a quantitative understanding of event-by-event fluctuations a much better understanding 
of the physics of the freeze-out stage is needed (see e.g. a recent discussion in Ref. ||2^). 

5 Summary and Final Remarks 

We have discussed event-by-event transverse momentum fluctuations in non-equilibrium 
and equilibrium models of nucleus-nucleus collisions in terms of the variable ^t- We have 
shown that the non-equilibrium initial state scattering models lead to a large (up to a 
factor of 10) enhancement of fluctuations for central Pb+Pb collisions at the CERN SPS in 
comparison to that observed in p-l-p interactions. 

The preliminary results of the NA49 Collaboration on fluctuation measure for central 
Pb-|-Pb collisions at 158 A GeV indicate a significant reduction of these fluctuations relative 
to fluctuations in p-l-p interactions . Thus the data show qualitatively different behaviour 
than that predicted in the initial state scattering models. We conclude therefore that such 
models need to be significantly modified in order to describe collisions of heavy nuclei. 

We have argued that even in the extreme case of global thermal equilibrium reached at 
the early stage of the collision the values of <I>t should be larger than zero (independent 
particle emission limit) due to the presence of thermal fluctuations at the early and freeze- 
out stages and the correlation between the multiplicity (entropy) and transverse momentum 
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(temperature) in the equilibrated massless gas. However as the origin of fluctuations is very 
different for equihbrium and non-equihbrium cases a significant change of the value is 
expected to occur in the course of an increasing equilibration of the created matter. 

We hope that the first exciting results of the NA49 Collaboration will prompt further 
extensive theoretical and experimental studies. In particular the measurements of event -by- 
event fluctuations in p+p, p+A interactions and central A+A collisions at various collision 
energies are necessary. 
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